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Abstract 

The purpose of this paper is to find explicit formulas for basic objects pertaining the 
local potential theory of the operator (/ — A) a / 2 , < a < 2. The potential theory of 
this operator is based on Bessel potentials J a = (I — A)~ a / 2 . We compute the harmonic 
measure of the half-space and write a concise form of the corresponding Green function 
for the operator (7 — A) a / 2 . To achieve this we analyze the so-called relativistic a-stable 
process on M. d space, killed when exiting the half-space. In terms of this process we are 
dealing here with the 1-potential theory or, equivalently, potential theory of Schrodinger 
operator based on the generator of the process with Kato's potential q = — 1. 
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1 Introduction 



The so-called interpolation spaces play an important role in harmonic analysis and partial 
differential equations (see, e.g. [5] and [H]). The most classical ones are Sobolev spaces and 
spaces of Bessel potentials. While (fractional) Sobolev spaces L p a are defined in terms of Riesz 
potentials I a = (— A)~ Q / 2 , the latter ones employ Bessel potentials defined as J a — (I — A) - "/ 2 . 
As Stein pointed out in his monograph [S], the both potentials exhibit the same local behaviour 
(as |x| — > 0) but the global one (as \x\ — > oo) of J a is much more regular. It is remarkable that 
both potentials can be analyzed in terms of stochastic processes (at least for < a < 2); the 
Riesz potentials are closely related to a-stable rotation invariant Levy process while the Bessel 
potentials, in turn, can be investigated in terms of the so-called relativistic process. 

Potential theory based on Riesz kernels (or, equivalently: potential theory for a-stable 
rotation invariant Levy process) is well developed and rich in explicit formulas, much like in 
the classical case of Brownian motion process. The homogeneity of Riesz kernels yields many 
elegant and transparent formulas for harmonic measure and Green function for such basic sets 
as balls and half-spaces in M. d (see e.g. |BGRj ). These formulas played an important role in 
setting up the so-called boundary potential theory of the operator (— A)"/ 2 and the Schrodinger 
operator based on it (see e.g. [BB1] and [BB2j or |CSlj ). 

In contrast to this situation, up to now, there were no explicit formulas known either for 
harmonic measure or Green function for the relativistic process for sets such as half-planes or 
balls. Nevertheless, an adequate boundary potential theory (at least for bounded smooth sets) 
was set up by Ryznar |Ry| . Let us point out that in recent years a number of publications 
concerning the potential theory of relativistic process and in particular Green function of this 
process appeared QCS2] . [K| . |KL] . |GRy| ); all however restricted to bounded sets. 

In this paper we provide the explicit formulas for harmonic measure and Green function for 
half-spaces for the operator — (/ — A)"/ 2 . Our approach consists of considering the operator 
H a = I — (I — A)"/ 2 , which is the infinitesimal generator of the relativistic a-stable process 
and regarding the operator — (/ — A)"/ 2 as the Schrodinger operator of the form S a = H a + ql, 
with q = —1. We then identify Bessel kernels as 1 -potentials for the relativistic a-stable 
process and compute the corresponding 1-Poisson kernel and 1-Green function for half-spaces. 
It is remarkable and surprising that even though our kernels do not exhibit either form of 
homogeneity, the resulting formulas are very transparent and very similar to these for a-stable 
case. 

The organization of the paper is as follows. We compute first the formulas for harmonic 
measure and Green function for one-dimensional case. To do this we apply complex-variables 
methods and some real-variables manipulations with definite integrals to obtain a satisfactory 
form of Green function. The d- dimensional case is then settled via application of (d — 1)- 
dimensional Fourier transform. For technical reasons, we have to consider the Poisson kernel 
and Green function not only for the operator — (J — A) Q / 2 but also for —{m 2 / a I — A) Q / 2 . 
Let us point out that we do not apply Kelvin's transform (which is the indispensable tool in 
the a-stable case). When performing a suitable limiting procedure we obtain the well-known 
formulas for the a-stable case. The last section is devoted to various estimates for the Green 
function of the half-space H, computed for the relativistic process (that is, for the operator 
H a = I — (I — A)"/ 2 ). To distinguish it from the corresponding object computed for the 
operator —(I — A) Q / 2 we call it O-Green function. The estimates are precise for x, y E M such 
that | a; — y\ < 1. To the best of our knowledge, it is the first case when the Poisson kernel and 
Green function of the relativistic a-stable process of an unbounded set are examined. 
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2 Preliminaries 



Throughout the paper by c,C,Ci ... we denote nonnegative constants which may depend on 
other constant parameters only. The value of c or C, C\ ... may change from line to line in a 
chain of estimates. 

The notion p(u) ~ q(u), u G A means that the ratio p(u)/q(u), u G A is bounded from 
below and above by positive constants which may depend on other constant parameters only. 

We present in this section some basic material regarding the a-stable relativistic process. 
For more detailed information, see Ry and [Cj. For questions regarding Markov and strong 



Markov property, semigroup properties, Schrodinger operators and basic potential theory, the 
reader is referred to [ChZ] and |BGj . 

We first introduce an appropriate class of subordinating processes. Let 9 a (t,u), u, t > 0, 
denote the density function of the strictly a/2-stable positive standard subordinator with the 
Laplace transform e~ txa/2 , < a < 2 . 

Now for m > we define another subordinating process T a (t, m) modifying the correspond- 
ing probability density function in the following way: 

6 a (t, u, m) = e mt 6 a (t, u) e- m2/au , u>0. 
We derive the Laplace transform of T a (t,m) as follows: 

E e -XT a (t,m) _ e mt ^(A+m 2 ^) / 2 _ ^ 

Let B t be the symmetric Brownian motion in IR d with the characteristic function of the form 

E o e i£-B t = e -t\e . (2 ) 

Assume that the processes T Q (t,m) and B(t) are stochastically independent. Then the process 
Xt' m = BT a (t,m) is called the a-stable relativistic process (with parameter m). In the sequel we 
use the generic notation X™ instead of J"'™. If m = 1 we write T a (t) instead of T a (t,m) and 
X t instead of X^. 

X™ is a Levy process (i.e. homogeneous, with independent increments). We always assume 
that sample paths of the process X t m are right-continuous and have left-hand limits ("cadlag"). 
Then X™ is Markov and has the strong Markov property under the so-called standard filtration. 

Various potential-theoretic objects in the theory of the process X™ are expressed in terms 
of modified Bessel functions K v of the second kind, called also Macdonald functions. For 
convenience of the reader we collect here basic informations about these functions. 

K v , v G M, the modified Bessel function of the second kind with index z/, is given by the 
following formula: 

r°° 2 

K u (r) = 2~ 1 ~ v r u / e~ v v^'" dv , r > 0. 
Jo 

For properties of K v we refer the reader to [E]. In the sequel we will use the asymptotic 
behaviour of K v : 

K v {r) = ^-T-Y r^0 + , u>0, (3) 



2 V2, 

K (r) = -logr, r^0 + , (4) 



K u (r) = ^=e~\ r^oo, (5) 
'It 
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where g(r) = f(r) denotes that the ratio of g and / tends to 1. For v < we have K v (r) = 
K-v{t), which determines the asymptotic behaviour for negative indices. 

The a-stable relativistic density (with parameter m) can now be computed in the following 
way: 

poo 

PT( X ) = / 6 a (t,u,m) g u (x)du, (6) 



where g u (x) is the Brownian semigroup, defined by ([2D- A particular case when a = 1 is called 
the relativistic Cauchy semigroup on M. d with parameter m. The formula below exhibits the 
explicit form of this density: 

Lemma 2.1 (relativistic Cauchy semigroup). The density p™ of the relativistic Cauchy process 
is of the form: 

P T(x) = 2(W2 7 r)^/ 2 ^^ +1)/2(m(|x|2 j t2)1/2) . 



(|x| 2 + t 2 ) 4 

Proof. We carry out the corresponding computations as follows: 



(?) 



OO 1 . 

1 I™|2 I a,, _™2„, i _Q/o _+2 / 



(47T U ) d / 2 y/A^ 

te mt f°° m 2„ /]„]2,4.2\ I*., du 



e -m 2 u e -(\x\ 2 +t 2 )/4u_ 



/j \±t± 4+1 I 1 

(4-7TJ 2 JO U 2 

2(W2vr)^)/ 2 ^^ +1)/2(m(|x|2 d + / )1/2) - 

(|x| 2 + t 2 )~ 



□ 



In the case of arbitrary a, < a < 2 we have the following useful estimate (see 
Lemma 2.2. There exists a constant c = c(a,d,m) such that 

maxp™(x) < c{t~ d ' 2 + r d ' a ) . (8) 



In the next lemma we compute the Fourier transform of the transition density (J6]): 

Lemma 2.3 (Fourier transform of p™). The Fourier transform of a- stable relativistic density 
p™is of the form: 

pfiz) = e mt e- t ^ 2+m2/a)a/2 . (9) 

Proof. 

p™(z) = [ p r t n (x)e i{z > x) dx = [ [ e mt g u (x)e- m2/au e a (t,u)due i{z ' x) dx 



R d JO 

OO POO 



e mt I e 



e mt e -t{\z\ 2 +m 2 / a ) a / 2 



□ 
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Specifying this to the case a = 1 we obtain 

From the form of the Fourier transform we have the following scaling property: 

p™(x) = m d / a P l rnt (w}/ a x). (11) 

In terms of one-dimensional distributions of the relativistic process (starting from the point 0) 
we obtain 

Xm 1/a vl 

where X™ denotes the relativistic a-stable process with parameter m and " ~" denotes equality 
of distributions. Because of this scaling property, we usually restrict our attention to the case 
when m — 1, if not specified otherwise. When m = 1 we omit the superscript "1", i.e. we write 
Pt(x) instead of p\{x). 

We work here within the framework of the so-called X-potential theory, for A > 0. The 
standard reference book on general potential theory is the monograph |BGj . For convenience of 
the reader we collect here the basic information with emphasis on what is known (and needed 
further on) about the a-stable relativistic process. 

We begin with the kernel of the resolvent semigroup of the process X™. We call the kernel 
corresponding to parameter A the X-potential of the process and denote it by U™(x) . It is of 
a particular simple form when A = m and we state it for further references. Again we denote 
by U\(x) the A-potential for the case m — 1. 

Lemma 2.4 (m-potential for relativistic process with parameter m). 

U™{x) = C(a, d) m - | g ( d-a)/2 > ( 12 ) 

where C{a,d) = r{a/2)nd/2 ■ 

Proof. We provide calculations for m = 1 and the general case follows from (1111) . 

poo /»oo /*oo 

Uiix) = I e~ t p t (x)dt = / / g u (x)e~ u 9 a (t,u)dudt 
Jo Jo Jo 

e-4r e - u [\ 6 a (t,u)dt)du 



(47ra) d/ 
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1 f°° _„ du 

e e 



r(a/2)(47r)^y 

2 l^ d+ a)/2 Kid _ a)/2 ) 

T(a/2)TT d / 2 \x\( d ~ a )/ 2 ' 



□ 



The first exit time of an (open) set D C M d by the process X™ is defined by the formula 

r D = mf{t > 0; X t m ^ D} . 

The basic object in potential theory of X" 1 is the X-harmonic measure of the set D. It is 
defined by the formula: 
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P£(x,A) = E x [td < oo;e-^l A (X™)}. (13) 

The density kernel of the measure P^(x,A) (if it exists) is called the X-Poisson kernel of the 
set D. 

In Section 3 and 4 we prove the existence of the m-Poisson kernel providing at the same 
time an explicit formula for it. The more general case of existence of A-Poisson kernel can be 
deduced from papers [IW] and [St] . 

Another fundamental object of potential theory is the killed process X^ ,D when exiting the 
set D. It is defined in terms of sample paths up to time td- More precisely, we have the 
following "change of variables" formula: 

E*f{X?' D ) =E x [t< t d - f(Xp)} ,t > . (14) 

The density function of transition probability of the process X t m ' D is denoted by p^ ,D . We have 

p?' D (x, y) = pT(x -y)-E x [t> r D ; pt TD (X™ -y)], x,yER d . 

Obviously, we obtain 

p?' D (x,y)<p?(x,y), x,yER d . 

p™' is a strongly contractive semigroup (under composition) and shares most of properties 
of the semigroup p™ . In particular, it is strongly Feller and symmetric: p^ ,D (x, y) = pT ,D (y, x). 
When m — 1, we write, as before, pf, instead of p\ L '. 

The A-potential of the process X t m,D is called the A- Green junction and is denoted by G^. 
Thus, we have 

rod 

G x D (x,y)= / e- xt p?> D (x,y)dt. 
Jo 

The " sweeping out" procedure provides another formula for the A- Green function of the set 
D in terms of the A-harmonic measure: 

G x D {x iV ) = Uf- f U?(z,y)P x (x,dz). (15) 

We now state some basic scaling properties both for the m-Poisson kernel and the m-Green 
function. The proof employs the scaling property ffTTj) and consists of elementary but tedious 
calculation and is omitted. 

Lemma 2.5 (Scaling Property). Let D be an open subset of M. d and PJj, G 1 ^ be m-Poisson 
kernel, or m-Green function, respectively, for D . Then 

P%(x, u) = m d l a P l ml/aD (m l / a x, m l ' a u), x^D,u^D c , 

G%(x, y) = m^/ a G l ml/aD {m l ' a x, m l ' a y), x E D,y E D . 
Thus, if D is a cone with vertex at we obtain: 

P£(x, u) = m d/a P l D {m l/a x, m l/a u), xED,uED c , 
GP£(x, y) = m^ a)/a G l D {m l/a x, m 1/a y), x E D,y E D . 
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Taking into account the above lemma and the fact that we restrict our attention to the case 
when D — H, the upper half-space, we usually assume that A = m = 1 and, as before, omit 
superscript "1", that is, we write Pd(x,u) and Gd(x,v) instead of Pp(x,u) and G\,{x, y). 

We also recall the form of the density function v m [x) of the Levy measure of the relativistic 
a-stable process (see 



Lemma 2.6 (Levy measure of relativistic process with parameter m). 



7r d / 2 r(l - |) V \x 



d+a 



» m ( x ) = ^Fn a\ ( \ — r~ ) K^m^lxl). (16) 



When m = 1 we omit, as usual, the superscript "1". 

The next lemma is taken from [G] . For reader's convenience we provide its proof. 
Lemma 2.7. For m = 1 we have 

P%m>t)<C^±M y t>2,x d >0. (17) 

Proof. Let Y t = X\ , where X t = (X t , . . . ,X^). By the symmetry of the random variable 
Y t we obtain 

P x (r m >t) = P x (infY s > 0) = P°(inf(-y a + x d ) > 0) = P°(supY s < x d ). 
Using a version of the Levy inequality ([E], Ch.7, 37.9) we have for any e,y > that 

n 

2P°(Y t ^y + 2e) - 2V P°(F« - Y Kk -p ^ e) ^ P°(sup Y* ^ y). 

k=0 Ki = n 

Note that X]fc=o-^ , °(^— — Y t (k-i) ^ e) — nP°(Y± ^ e) — > t j°° is(x)dx, hence, by symmetry 

71 n TL 

again 

poo poo 

P°(\Y t \ ^y + 2e)-2t / u(x)dx = 2P°(Y t ^ y + 2e) - 2t / ^(x)rfx ^ P°(sup y ^ y). 
This implies that 

POO 

P x ( m >t) = P°(supy < x d ) < P°(\Y t \ <x d + 2e) + 2t / 
For £ ^ 1 we obtain from (1161) and §5§ 

v{x)dx ^ Ce- £ e~ a/2 . 

Because of the Lemma 12.21 we have that the density of Y{t) is bounded by Ct" 1 / 2 , t>2 hence 
taking e = | hit we obtain 

P x {m>t) ^C(x d + in t)r 1 ' 2 . 

□ 

In order to improve the above estimate for x close to the boundary we use the following 
result proved recently in |GRy| . 
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Lemma 2.8. Assume that d=l. Let D = (0, 1) and x G D. Then 

(x(l - x)) a/2 

P X (X TD G dz) « ^| r e~\ x E D, z > 1. 

[z — l) a ' z (z — X) 

This implies that 

E X [X TD > 1;X TD ] « P X (X TD > 1) « x Q / 2 . 

VKe a/so /lave t/iat 

E x td w (x(l -a;)) Q/2 . 

Actually in |GRy| it was shown that the Green function of D is comparable with the Green 
function of the corresponding stable process. By standard arguments (see |Ry| ) this implies 
the above lemma. 

We further need the following strenghtening of Lemma 12.71 

Lemma 2.9. For < Xd < 1/2 we have 

P x (m >t)< Cx a J 2 lnt/t 1/2 , t > 2, (18) 

where C is a constant. 

Proof. It is enough to prove the claim for d = 1. Let D = (0, 1) and assume that < x < 1/2. 
By the Markov property and then by Lemma [2.71 we obtain for t > 2: 



^'(m > t) = P x (t d >t,r D = m) + £*[t d < r e ; (r H > *] 

< P x (rD >t) + E x [t d < m; X TD + \nt]/t 1/2 

= P x (r D >t) + E X [X TD > 1 ; X TD \lt x l 2 + hit P X (X TD > l)/t 1 ' 2 

< Cx^lnt/t 1 / 2 . 

The last inequality follows from Lemma 12.81 The proof is complete. □ 



3 Poisson kernel and Green function of — (m 2//a 7 — j^) a ^ 2 
for M+ 

We first consider the one-dimensional case. For technical reasons (see Section 4) we have to 
distinguish notation in this case from the d-dimensional one so we denote the relativistic Poisson 
kernel by Q 1 ^^ (instead of P|f as previously). 

Theorem 3.1 (one-dimensional relativistic Poisson kernel). Denote 

E x [e- m ^-X™ Q x) E du] = Q? 0tOo) (x,u) . 

Then we have for u < < x: 

y ' ' 7T \ — U I X — U 
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Proof. By the scaling property it is enough to deal with m — 1. According to the general facts 
of potential theory (see e.g. |BGRj ) it is enough to verify that the measure Q(o,oo)(^, u) satisfies 
the following "sweeping out" principle: 

<3(o,oo)0, u)Ui(u, y) du = Ui(x, y) , 

for all y < 0, where Ui(x,y) is the 1-potential of the relativistic a-stable process (see ( fl2|) ). 
This is, however, shown in the following lemma below. □ 

Lemma 3.2. For x > > y we have 



sin 



o 



x 

-u 



-\x-u\ Ki- a (\u — y\) Ki-a (\x — y\) 



i 1 1 ~ a 
\u — y\ 2 



du 



i 1 1 ~ a 
\x — y\ 2 



7T J \—U J \X — U\ 

Proof. Let x > > y and consider the following function of complex variable z: 

1 e z ~ x Ki^*(z - y) 



l-a 
2 



Z 2 x-z (z-yY 

This function is holomorphic in C\(— oo, 0]\{x}. We are going to integrate the function above 
over a contour described below. We make the branch cut along the axis (— oo, 0] and make the 
contour of integration to wrap around this line (see picture) and we add a circle around x say 

73- 



-»■ 71 


-Q 1 


y 72 


73 



By the Cauchy theorem, we get 



Using the asymptotic expansion for the modified Bessel function K v (z) we get 



(19) 



Kx=<* (z - y) 



( z ~y) 



l-a 
2 



(z-y) 

9 



— ^R (z - y), 



where Rq(z) = 0(1) and z is large enough. Using the expression given above it is easy to show 
that the integral over T of f(z) vanishes when r — > oo. 

To calculate the integrals over 71 and 72 we examine the behaviour of the function / near 
the branch cut. We have the following relations for the modified Bessel functions for u < y 

lim Ki-a (u — y + ie) = ~Ki- a (\u — y\) — inli-a (\u — y\), 

£ — ^O - ! - 2 2 2 

lim Ki- a (u — y — ie) = Ki-*(\u — y\) + mlx- a (\u — y\). 

t— >o+ 2 22 



We have also 



lim (u — y + ie) 
lim (u — y — ie) 



1-q 
2 



1-q 
2 



Z7r(q — 1) _ 1 — q 

e 2 |u — y| 2 . 

i7r(l — q) _ 1 — q 

e 2 |u — y\ 2 , 



Next, for w < we have 



lim (u + ie) 2 

e-»0+ ' 

lim (u — ie)~ 2" 



e 2 kt 2 . 



e 2 hi 2 



Using all the relations given above we obtain that for u < y 

A"l_q {\u~y\) Ii- a (\u-y\) 



fi( u ) = lim e-*o+ f(u + ie) 
/f(u) = lim e _+ + f(u - ie) 
Similarly for y < u < we get 



1 e"" 



|w|? x-u 
1 e"" 1 



? 7T CI 

e 2 



|u|7 



e 2 



1-q "I - 71" 1-q 

\u-y\~5~~ \u~y\~^~ 

Xi_q(|n-j/|) 7i_q(|n-y|) 
7T- 



W-y\~ 



/ 2 +(u) = lim e ^ + /(« + ie) 
/if (u) = lim e ^o+ f(u - ie) 



-e 2 



1 -U-3! i7ra _a (1^ 2/|) 

-e 2 =az . 



l"-2/l" 



Combining all above we find that 

I e -\x-u\ Ki-a (\u — y\) 



net \ />0 



7T 



-u)^ \x-u\ \u-y 



1-q 
2 



du. 



We have also 



2ttz 



77- / 



73 



1 Ki_^(\x-y\) 

ST ! 1 1-q ' 

x 2 \x-y\~ 



Using (Tl9l) and the relations given above we obtain the desired formula. 



□ 
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Theorem 3.3 (one-dimensional relativistic a-stable Green function). For x,y > we have 



GTo^y) = 1^1^ e -^/. l ,_ y|( , +1) v 2tQ/2 „ 1(t + irl/2 rft _ 

Proof. By the scaling property we may assume m = 1. Due to the symmetry of the Green 
function, it is enough to determine G(o i00 )(i, y) for < x < y. We compute the compensator of 
the Green function for the one-dimensional a-stable relativistic process. We want to evaluate 
the following expression 

H(x,y) = E'le-^-WdXr^yV)] 

= CM ^i f (±.y /2 €^i^z$ dUi (20) 

V J-Oo\-Uj X-U — 

where C(a,l) = w~j^r/2 • Substituting (— u) l ~ a l 2 = v and taking into account the following 
well-known identities 



v 2/(2-a) + x 



OO 







Ki-a(z) 1 r°° e ~ zt 



C((X,1) zk* ~ r(a/2)r(l - a/2) J, (t* - 1)^ dt 



we obtain that the right-hand side of ( 12011 is of the form 



(2 - a)7rr(a/2)r(l - a/2) y U J 7i (* 2 - l) a/2 

2sin(vra/2)x-/ 2 e- /"* /"» f°° , w+t+1)v 2 K 2- a) , e"* 



wx / e _ (w+t+1)v 2/(,- a) dv e ^ ^ 



(2 - a)7rr(a/2)r(l - a/2) J 1 J J (P - 

The interior integral can be expressed as 

OO POO 

JO 

= ^_^r(l-a/2)( U ; + t+l)( Q - 2 )/ 2 . 

Consequently, we get 

Observe that the expression in brackets can be computed as follows 

POO POO 

x a/2 e -x(v,+t+l)( w + t+1 ya-2)/2 dw = / 

JO Jx(t+1) 

s(t+l) 



T(a/2)- / e-'s^^ds 
Jo 



F{a/2) - ~x a/2 (t + 1) Q/2 / e" 1 " /ax(t+1 W 
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sin(vra/2) f°° e-fo-*)* ^ 2 sm(Tra/2)x a / 2 f°° f 1 e -(»-*)* 



Thus we get 

K ' Ji (t 2 -l)«/ 2 a7rT(a/2) J, J (t - 1)«/ 2 

Because 

G( ,oo)(x,y) = Z7i(z,y) -^[e-^-)^(X r(0oo) ,y)], 

thus we finally have 



2x a/2 


ari 


:«/2) 2 r(i - 


a/2) 




2x Q / 2 




aTi 


:«/2) 2 r(i - 


a/2) 




2x OL l 2 e x ~ y 




aTi 


:«/2) 2 r(i - 


a/2) 


2x Q 




-2xw 2 / 



(t - 

it - 1)<*/ 2 



1 /*oc „—(y—x)t 

e~ xw2/a / e -^ 2/Q ^ — dirfu, 



aT(a/2) 2 J 



e- 2xw2/a {xw 2 ' a + y-x) a ' 2 - l dw 



X e -2v v a/2-l^ + y _ x ^/2-l dv 

4xy 

e - iu+ (y-x)^ u */2-l {u +{y _ x)2) -l/2 ^ 



r(«/2)2 y 
1 

2«r(a/2)* 7 
2«r(a/2) 2 y 



□ 



4 Poisson kernel and Green function of —(m 2 l a I — A)"/ 2 
for M 

Our proof relies on the computations of the (d — l)-dimensional Fourier transform of the d- 
dimensional Poisson kernel as well as the corresponding d-dimensional Green function. Namely 
we show that these Fourier transforms can be expressed in terms of the corresponding one- 
dimensional objects, which we can easily invert. 

To avoid confusion, we introduce the following notation, to distinguish one-dimensional and 
c?-dimensional objects. 

Notation: 

xGl d_1 , x = (x, Xd) £ M. d . 
In the proofs below one-dimensional quantities play an important role, hence we denote by 

q™(x) — one-dimensional a-stable relativistic density with parameter m, 

V™(x) — corresponding A-potential, 

Q(0 )OO ) ( x ? u ) ~ corresponding Poisson kernel. 
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We begin with computing of the (d — l)-dimensional Fourier transform of the A-potential 
£/\((x, Xd), (y,Vd)) = U\((x — y), (x d — yd)) with respect to the variable y. To avoid nota- 
tional complexity, we denote it by U\(x, -){z). We employ this kind of notation throughout the 
whole section. Thus, we have 

Lemma 4.1. 

(Mx7-)(z) = e^Vf (xd-Vd), 
where k = (|z| 2 + iy/ 2 and A = n a + A — 1. Specifying to the case A = 1: 

lMx7-)(z) = e^V^ixd-Vd) 

-e Hz ' x) K 2 2 - ^ — . (21) 

\xd-Vd\ 2 



v^r(a/2) 



Proof. We begin with computation of the (d— l)-dimensional Fourier transform of the transition 
density of the normal distribution: 

g^)(z) = [ g u (x - y)e^dy = ^^e^Sr 1 ^. 
J W d-i {Anuy/ 1 

In the next step we use this to find the (d — l)-dimensional Fourier transform of the transition 
density p t : 



Pt{x,-){z) = / Pt {x-y)e i{z ' y) dy = e t g u (x)(z)e u 9 a (t,u)du 

jRd-l JO 



= ei (^) e *(i-(H»+i)-/») / ^L^ e t(|z| 2+1 )^ e -(| z | 2+ i) Ue -toL^ (t )du 

J0 (47Tu) 1 /2 

= e i(-.*) e -(*-i)t [°° A^-u 1 e _ (a^i ga(t , d 

Jo {A-ku) 1 / 1 

where m = (|z| 2 + 1) Q / 2 . Note that the the integral expression in the last line is the one- 
dimensional a-stable relativistic density with parameter m which we denote by q™. Hence we 
obtain 

p^v)(z) = e^ed"*^* (x d - y d ). 
As a consequence we obtain the (d — l)-dimensional Fourier transform of the A-potential : 



U x (x,-){z) = e^' x M e£-*- x \t(x d -y d )dt 

Jo 

= e^V&^ixd - y d ) = e^V™(x d - y d ), 

where A = (m + A — 1) and V~™ is the A potential for the one-dimensional relativistic semigroup 
q™ with parameter m. Now we take A = 1 then A = (rh + A — 1) = (|z| 2 + l) a / 2 = fh. Denote 
K = k(z) = (|z| 2 + l) 1 ' 2 . By (TjED we have 



x/Th/ \ T ««, v 2^ t^Kx^iKlxd-Vdl) 

V~™{x d - y d ) = V£ (x d - y d ) = k * 

y/TTl (Oi/l) 



_2 



□ 
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Our next step is to find the (d— l)-dimensional Fourier transform of the Poisson kernel and 
the Green function. 
Recall that 

P m (x, u) = E x [e~ m , X m edu], x6l,neff, 

is the Poisson kernel for the operator —(I — A)"/ 2 on the set EI. We also recall the formula (|15p 
for the Green function, applied for m — 1: 

Gm(x,y) = U 1 (x-y)-E x [e-^U 1 (X m ,y)}, 

for i,i/6E The function Gu{x, y) makes sense also for x G EI and i/Gff and from the general 
theory it is zero for this range of x and y. This actually gives a condition for the Poisson kernel: 

= G a (x,y) = U 1 {x-y)-E x [e^U 1 (X m ,y)} 

= U^x-y)- I P m {x,u)U 1 {u,y)du, (22) 



for x G HI and y G H c . Now we take the Fourier transform of Gu(x,y). We carry out com- 
putation for x = (0, x d ) since the general case can be easily deduced by translation invariance 
of the process. Denote k = k(z) = (|z| 2 + l) 1 / 2 and recall that V™ is the m-potential for the 
one-dimensional relativistic semigroup with parameter m. Then, 

G^TO(z) = U^-)(z)-E x ie-^U 1 ~(X^,-)(z)} 

= V£{x d - y d ) - E*[e—e^^V::{X d m - y d )}. 

Next, 

E x [e-^^) V ^{X d m -y d )} = [ e^P m {x,u)Vi::{u d -y d )du 

Ju d <0 

= [ f[ e l{z ' u) P m ((0,x d ),(u,u d ))du)v::(u d -y d )du d 

J -co VjR d ~ 1 ) 

r° - — - 

= / P m {x,-){z)V^{u d -y d )du d . (23) 

J — oo 

Hence the condition ([22]) is equivalent to: 

V*™(x d -yd)= / Pm(x, -)( z )V^(u d - Vd)du d , 

J —oo 

for y d < < x d . Now this implies by the one-dimensional result ( Theorem 13.11) that 



P H (x, -)(z) = Qto,oo)(^,u d ) = CI - f . (24) 

v ; \~ u dJ Fd - u d \ 

Note that we used the fact that Pm(x, -)( z ) is real by symmetry of the process. 

Next we proceed to computation of the Fourier transform for Gm(x, y) for x = (0, x d ),y G EL 
It immediately follows from ( 1231) . ( J24|) and Theorem 13.31 that 

G m {x, ■) (z) = ( V£* {x d -yd) - / Qm ( x d, u d )V^ (u d - y d )du d 



G(o,tx>)( x d, y d ) 



pAXdUd — 1 

-(«+(xd-yd) a ) 1/a (l«l a +i) 1/a ds. 



2«r(a/2) 2 7 (s + (x d -y d ) 2 ) 1/2 
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Theorem 4.2 (Poisson kernel). Let 

E*[e- m ™,X™edu] = P^(x,u). 

Then we have 

m i/<*Y /2 ( x d \ a/2 K d/2 {m l l a \x-u\) 



iff(*,«)=2C^— j (-i-j ^— ± u d <0< Xd . (25) 

Proof. We prove only the case m = 1 since the general case follows from the scaling property. 
Also it is enough to consider x = (0,x d ). Applying ( 1241) we have 



P H (Z, 0( z ) = Q(0,oo)( x d,U d 

= cl 



x \ Q / 2 g-K|a; d -n d | 



-u d J \x d -u d \ 
i f x rf V /2 e-d z l 2+1 ) 1/2 l^-^l 
a V-^d / \ x d - Ud\ 

Taking into account (fTUl) and ((7|) with (d— 1) instead of d and \x d — y d \ instead oft we complete 
the proof. □ 



2 l-a m d/2a\ x _ \a-d/2 



Theorem 4.3 (Green function for 

Proof. We prove only the case m — 1 since the general case follows from the scaling property. 
Also it is enough to consider x = (0,x d ). Recall that 

GH(X ' m ~ 2«I>/2) 2 7 {s + (x d -y d yy/> e dS ' 

Taking into account fTTUl) and ((7|) with (d — 1) instead of d and (s + (x^ — yd) 2 ) 1 ^ 2 instead of t 
we obtain for d > 1 

^ , v f AxdVd H- 1 ^/ 2 ((|y| 2 + (a: rf -l/ rf ) 2 + g ) 1 / 2 ) , 



(27r) rf / 2 I>/2) 2 y (|y| 2 + (x d - y d f + S )V 2 )<*/ 

S 2 -r, 7^ TT71 dS 



(2yr) rf / 2 r(a/2) 2 7 (\x - y\ 2 + s) d l A 



j^ 2 {t + ~iy/* Kd/2ilx ~ yl{t + 1)V2) dt ■ 



(27r) d / 2 r(a/2) 2 7 (t + l) d / 4 

Corollary 4.4. Assume that \x — y\ < 1. T/ien there is C = C(a,d) such that 

(x d y d Al) a/2 <CG m (x,y). 
This estimate is not optimal but sufficient for our purposes in the next section. 



□ 
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Proof. First, observe that K d / 2 (r)/r d ^ 2 is decreasing (see, e.g. [EJ). Next, we use one of the 
form of the Green function obtained in the above proof to arrive at 

G M [x,y) = , / — ds 



(27r) d / 2 r(a/2) 2 J (\x - y\ 2 + s) d l 

2^ [ 1M ™« q_ 1 ^ d/2 ((|x-y| 2 + 8 ) 1 / 2 ) 

" (2vr) rf / 2 r(a/2) 2 y 52 (|x - i/| 2 + s) d / 4 * 

2 1 -" /- 1A4 ^ h_ 1 A> 2 (2) 

> 7Z 777777 / S 2 TJIi — «s 



(27r) d / 2 r(a/2) 2 J 2 d / 4 
C{Ax d y d M) a/2 



□ 



Note that using the asymptotic behaviour of Kd/2{f) as r — > + we obtain the well-known 
formulas for the Green function and the Poisson kernel for the standard symmetric (rotation 
invariant) a-stable process: 



„ - x dVd a 1 

I \x~y\' 2 t~2 

lim G m (x, y) = C{a, d) / rft, x d , y d > 0. (27) 



Similarly, 



, a/2 



lim Pg{x, u) = C(a, d) [ — ) , rr, n d < < x d . (28) 

m->0+ \—Ud J \X — 111 



As far as we know, it is the first alternative proof of the Poisson kernel formula for the standard 
symmetric (rotation invariant) d- dimensional a-stable process without application of Kelvin's 
transform. 

As a corollary we compute now the formula for E z e~ m . Its importance, among other things, 
is due to the fact that it is harmonic for the operator (/ — A)"/ 2 . In probabilistic terms, it is 
harmonic for the Schrodinger operator based on the generator of our relativistic process with 
the potential q — — 1. 

Corollary 4.5. For any < a < 2 we have 

1 r°° 

E z e -m t a/2-X-t dt 

r(a/2) J Zd 

Proof. The proof consists of computing the mass of the Poisson kernel. It is obvious that we 
may assume that d — 1. 

2 

Substituting (—it) = and taking into account the following well-known identity 

2 



2 

v 2 ~ a + x Jo 



1(3 



we have, after changing order of integration: 



VT / V — U / Z — U 



JO 



r(a/2) 

This clearly ends the proof. 



z a ' 2 / (w + ir^e-^+^dw 
Jo 

9 r°° 

t a - l e~ t2 dt 



r(a/2) 7^ 
1 



□ 



5 Estimates of Green function of / - (/ - A)"/ 2 for M 

In this section we work under the assumption that m = 1. In this section we deal with the 
Green function for the set HI, corresponding to the operator I— (I — A) Q / 2 . Equivalently, it is the 
usual Green function (that is, 0-Green function) for the relativistic process X t . To distinguish 
it from the previously considered Green function for the operator —{I — A) Q//2 we denote it by 
Gr jj. Our objective is to establish some estimates for G^. The estimates will be sharp if x, y 
are close enough. We apply for this purpose the estimate for Gh contained in Corollary 14.41 
The following lemma will be useful in the sequel. 

Lemma 5.1. There is C such that: 

Pt{x -y)- Pt (x - y*) < pf(x, y) < C{r d ' 2 + r d ' a )P x {m > t/3)P y (r H > i/3), 
where y* = (y u . . . ,y d -i, -yd), x,y e EL 

Proof. We start with the upper bound. Since pf- (x, y) is a density of a semigroup and pf(x, y) < 
m&x zeR dp t (z) then we have 

$t(x,y)= / pf(x J z)pf(z J y)dz <maxp t (z) / pf (x, z)dz = maxp t (z)P x (r m > t). 

Next we repeat that argument to have 

pl(x,z)pf(z,y)dz <m a xp t (z)P x (m>t) [ pf(z,y)dz 
= maxp t (z)P x (T M > t)P y (r M > t), 

zSR d 

which proves the upper bound since max 2gK d pt(z) < C(tr d l 2 + t~ d / a ) (see Lemma E2D. 



pl(x,y) 
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To get the lower bound we use the subordination of the process to the Brownian motion: 
X t = B Ta( ty Then 

pf(x, y) = P x (B Ta{t) e dy, B Ta(s) e H, < s < t) > P x (B Ta{t) e dy, B s e H, < s < T a (t)) 

Using the independence of T a and B we obtain 

P x (B Ta[t) edy,B s eU,0<s< T a (t)\T a (-)) = g Ta{t) (x-y)-g Ta{t) (x-y*), 

by the well known formula for the density of the killed Brownian motion on exiting H. Inte- 
grating and using the fact that E°gT a (t){z) = Pt(z) we obtain the lower bound. □ 

Lemma 5.2. There is C such that: 

G° m (x,y)>CG 9 m (x,y), 
where G^(x,y) is the Green function o/H for the Brownian motion. 

Proof. Let V(x, y) = f^°(pt(x — y) — p t (x — y*))dt. From the previous lemma it is enough to 
prove that V(x, y) > CG^(x,y). We have 



V(x,y) = / (pt(x - y) - jk(x - y*))dt 





oo POO 

e* / {g u {x -y)- g u {x - y*))e~ u 8 a (t, u)dudt 
Jo 

3 POO 

(g u (x - y) - g u (x - y*))e~ u I e t 6 a (t,u)dtdu 

Jo 

(g u (x -y)- g u (x - y*))G(u)du, 

where G{u) = e~ u J °° e t 9 a (t,u)dt is the potential of the subordinator T a {t). It was proved 
in |RSVj that G{u) is a completely monotone (hence decreasing) function and m.i u>0 G{u) = 
lim^oo G(u) = C. We find the constant C = lim^oo G(u) by taking into account the asymp- 
totics of the Laplace transform of G(u) at the origin: 

OO POO POO POO -| r) 

e~ Xu G(u)du= e* e~ u{1+x) 9 a (t,u)dudt = e'e'^+^dt 



JO Jo Jo 

Applying the Tauberian Theorem we obtain that C = 2/ a. 
Thus, since g u (x — y) — g u (x — y*) > we finally obtain 



(1 + A) Q / 2 - 1 Xa 



V(x,y) = / (g u (x - y) - g u (x - y*))G(u)du 
Jo 

POO 

> C / (g u (x - y) -g u (x - y*))du = CG 9 M (x,y). 



o 



□ 



The next result provides a general upper bound for the Green function. 
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Theorem 5.3. There is a constant C such that 

G° m (x, y) < C[<iP(x d )iP(y d ) + G H (x, y)], x, y e H, (29) 
where the function ip has the following form depending on the dimension d: 



for d = 1 
ford = 2 
for d > 3 



ip(y) = v a l 2 , < v < 1 and ip(v) = v 1 / 2 , v > 1 ; 
ijj( v )= v a / 2 i < u < 1 and ip(v) = ln 1/2 (1 + v) , u > 1 ; 
^) = («Al) a / 2 , w>0. 



Proof. The proof will rely on estimates of P x (tm > t) and the application of Lemma [5.1[ 
We proceed to estimate the Green function. First we split the integration 

> i*2 poo 

pf(x,y)dt = / pf(x,y)dt+ / pf(x,y)dt. 

JO J2 

To estimate the first integral we use 

! ,2 

pf(x,y)dt<e 2 / e- t pf(x,y)dt<e 2 G m (x,y). 
Jo Jo 

Due to Lemma [5. II 

POO POO 7, 

y tf(x iy )dt< J P x (m>t)py( m >t)^ = R(x,y). 

For d > 3 we have 



fl(x,y) < P*(r e > 2)P^(r H > 2) jf ^ < C(*„ A l)"' 2 (y d A 1) Q / 2 , 

due to ( fl8l) . which completes the proof in this case. 

To deal with d = 1, 2 observe that by the Schwarz inequality 

R 2 (x,y) < R(x,x)R(y,y), 

so we estimate R(x,x). First consider the case < a/2- Then using (ITTjl and (JT51) we have 

= I ( PX (tm > t)f W < Cx« ^ K —^ Wr 
If x d > a/2, using (fT7|) we estimate 

(^>t)) 2 ^ = / (P x (r m >t)f W2 +l (P%m>t)f W2 



.1 



t d / 2 J r 2 I t 1 / 2 J t d ' 2 ' 

Hence for d — 1 we have R(x,x) < Cx, while for d = 2we arrive at R(x,x) < Clna^. Taking 
into account all cases we get 



pf(x,y)dt<C^(x)^(y). 
The proof of the theorem is complete. 



□ 
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Theorem 5.4. For d = 1 and \x — y\ < 1, we have that 

G° m (x,y)~G m (x,y) + xAy. (30) 
For d = 2 and \x — y\ < 1 we have that 

G° m (x, y) « G H (x, y) + ln(l V (x 2 Ay 2 )). (31) 

For d > 3 we have that for \x — y\ < 1, 

G° m (x,y)^G m (x,y). (32) 

Proof. Step 1. We first show that if \x — y\ < 1 then 

G H (x, y) w Ge(x, y) if x d Ay d ^2 or if d ^ 3. (33) 
Recall that by Corollary 14.41 we have 

(x d y d Al) a/2 <CG m (x,y). (34) 

We first consider the case x d A y d < 2. The condition \x — y\ < 1 implies that x d V y d < 3. 
Hence 

XdVd < Q(lAx d y d ). 
Then due to Theorem 15.31 and (|3"4"|) we obtain 

Gm(x,y) < G° B (x,y) < C x {G m (x,y) + (x d y d ) a/2 ) < C 2 G m (x,y) . (35) 

We now consider the case when d > 3. It is elementary that 

(x d Al)(y d Al) < x d y d A 1 . 

Then the inequality (1351) can be rewritten in the following way regardless on the assumption 

on x d A y d , 

G m (x,y) < G° m (x,y)<C(G m (x,y) + ((x d Al)(y d Al)r/ 2 ) 

< C(G u (x,y) + (x d y d Al)^ 2 ) 

< C , iG H (x,?/), 

where in the last line we applied again (1341) . This completes the proof of Step 1 and proves the 
theorem for d > 3. 

Step 2. In this step we complete the proof of the case d — 1. 

The lower bound follows from the estimate proved in Lemma [5\2l : G° ^(x, y) > CG 9 ^ ^(x, y), 
where G 9 ^^ is the Green function of (0, oo) for the Brownian motion. 
It is well known that G 9 ^ (x,y) = x Ay so 

x A y + G (0 ,oo) (x, y) ^ (7G[ oo) (x, y) , 

and the lower bound is established. In the case x A y < 2 we also get the upper bound since 

G °(o,oo)( x ^y) ~ G (o,oo)(a;,y)) by Step 1. 
If x A y > 2 we obtain by Theorem 15.31 

G° { o,oo)(^y) < C 2 (G i0iOo) (x : y) + (xy) 1 / 2 ). (36) 
20 



Since \x — y\ < 1 and x A y ^ 2 we have (xy) 1 ^ 2 ^ x \/ y < (x A y) + 1 ^ 2(x A y), which 
completes the proof of ( l30l) for the case x Ay ^ 2. 

Step 3. Now we deal with the case d = 2. We claim that G^(x, y) ^ Cln(l V (a?2 A 7/2))- It 
is enough to show it for x 2 A y 2 > 1. It is well known that G^(x, y) — ^ In t^t ■ If |^ — y\ < 1 
then 

]_ I /jj* — y| ]_ 

G h(^*/) = 7T ln 1 r>^-^\x*-y\ 

2ir \x — y\ 2tt 

> ^-\n(x 2 + y 2 ) > — ln(l V (x 2 Ay 2 ))- 
As a consequence of the above inequality and Lemma 15.21 we obtain 

G m (x, y) < G m (x, y) + ln(l V (x 2 A y 2 )) < C(G m (x, y) + G 9 m {x, y)) < G° (2, y) (37) 

If now x 2 A y 2 < 2 then by Step 1 G^(x, y) ~ G^(x, y) and we obtain the conclusion. 
We proceed with the case x 2 A y 2 > 2. Note that x 2 V y 2 < 1 + x 2 A y 2 , so 

ln 1/2 (l + y 2 ) ln 1/2 (l + x 2 ) < 2 ln(x 2 A y 2 ) = 2 ln(l V (x 2 A y 2 )) 

Hence by Theorem 15.31 

G° m (x, y) < C(G m {x, y) + ln 1 / 2 ^ + y 2 ) hi l ' 2 {\ + x 2 )) < C(G m (x, y) + ln(l V (x 2 A y 2 ))), 

which, together with the estimate 0371) completes the proof of Step 3 and of the theorem. 

□ 

Remarks. 

1. The above theorem shows that, in particular, for fixed y G EI the function G^(x,y) ~ 
x a J 2 , for x near the boundary. With an extra effort we are able to strengthen this result and 
actually show that for fixed y G EI we have 

lim G° u (x,y)/xf = C(y,0. 

2. We are able to show that for d = 1 the estimate for G^ stated in Theorem 15.41 is optimal 
for the whole range of x, y G EI, that is, we have 

gU x iv) ~ G ^( x iy) + x A v 1 x ,y em. 
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